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The dispersion relation of longitudinal electrostatic oscillations in a relativistic plasma is studied in 
^ \ the context of the nonextensive statistics formalism proposed by Tsallis [C. Tsallis, J. Stat. Phys. 52, 

f">) ■ 479 (1988)], where nonextensivity is characterized by a parameter q in Tsallis's entropy, q — 1 

■ corresponds to the usual Boltzmann-Gibbs, extensive statistics formalism. In the nonrelativistic 

regime, normalizability of the equilibrium distribution function implies that — 1 < q < do. We 
^ ■ show that in the relativistic regime much tighter constraints must be satisfied, namely < q < 

1 + fesT/mc 2 , where ks is the Boltzmann constant, T is the temperature of the plasma, and m is 
the particle mass. 

C^l , Then we study longitudinal oscillations in a proton-electron plasma, assuming immobile protons, 

and electrons whose distribution function maximizes Tsallis's entropy. The dispersion relation of 
these oscillations is written in integral form for the long wavelength limit. Explicit expressions 
in terms of generalized hypergeometric functions can be found for all possibles values of q in the 
ultra-relativistic regime. 
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£3 ■ I. INTRODUCTION 

a 

f Traditionally, the equilibrium of statistical systems has been studied based on the Boltzmann-Gibbs entropy, 

Sbg = -kB } Vi In Pi , (1) 

>V 

■ where ks is the Boltzmann constant, and pi is the probability of the i-th microstate. However, an increasing amount 
|' of experimental, computational and theoretical evidence, shows that this formalism is not adequate for systems where 

long range interactions in time or space are important, and, in general, systems where the relevant phase space has 
fractal or multifractal structure. Tsallis [ij proposed the following form for the entropy of such systems: 

g: s^ks 1 - ^ , (2) 

■ where q is a real number. For q = 1, S q reduces to Sbg- 

S q has the basic property that for independent systems A and B, the entropy of the composite system A + B is 
^ , given by: 

S q {A + B) = S q (A) + S q (B) + (1 - q)S q (A)S q (B) . (3) 

Thus, q is a measure of the nonextensivity. Among the various physical systems where connections with the Tsallis 
^~>| entropy formalism have been found are gravitational systems 0, long range Hamiltonian systems Q, nonlinear 

■ inverse bremsstrahlung absorption in plasmas turbulence 0, and granular systems 0, and many applications in 
astrophysics, fluid dynamics, biology, economy, quantum mechanics, and others. (See, e.g. 0-) 

An important property of the nonextensive formalism is that the distribution function which maximizes S q is non- 
Maxwellian HQ. Specifically, for q < 1, high energy states are more probable than in the extensive case; for q > 1 
high energy states are less probable than in the extensive case, and there is a cutoff beyond which no states exist. For 



X 



a unidimensional gas where E = mv 2 /2, this cutoff is given by 



IW = y/2k B T/m(q - 1) . (4) 

Velocity distributions in plasmas are often found to be non-Maxwellian. High energy tails in space and laboratory 
plasmas |ld llll ll'2t ll.l |14| are a common occurrence. Metaequilibrium non-Maxwellian radial density profiles in 
cylindrical turbulent pure electron plasmas are observed [l5|. As stated above, such distribution functions emerge 
naturally if the system is described by a nonextensive statistics, and several authors have successfully applied this 
approach to plasma systems 0, [13, 13 ■ 

Lima et al. studied the longitudinal oscillations of a collisionless thermal plasma, in the context of the nonexten- 
sive formalism based on They derived the dispersion relation for electrostatic waves, and calculated it explicitly 
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for an equilibrium distribution function which maximizes • They found that it fits experimental data better than 
if the energy distribution is given by a Maxwcllian. 

All the previous analysis have been made for non-relativistic plasmas. However, relativistic extensions of the Tsallis 
formalism are possible |l9j. Our purpose is to study the consequences of a relativistic nonextensive statistics in 
plasmas. In particular, we will examine some general properties of the equilibrium distribution function and calculate 
the dispersion relation for electrostatic waves in the long wavelength limit for ultrarelativistic plasmas. 



II. RELATIVISTIC DISTRIBUTION FUNCTION 

In the nonextensive description, the equilibrium distribution function for a relativistic plasma can be written ^ 



Up) = A q 



i-(?-i) 



E 



k B T 



1/(9-1) 



where 



E = y m 2 c 4 + p 2 c 2 , 



(5) 



(6) 



is the particle energy and A q is a normalization constant. 

Equation © generalizes the nonrelativistic expression (see, e.g., ||). For q 
distribution function f q= i — Aiexp^—E/ksT). 

For simplicity, in the following we will consider a one dimensional plasma, so that p — pz. 

Since the relativistic energy is positive, f q will be real and normalizable if E satisfies the following conditions: 



if q > 1 
if q < 1 



< E < 



0<E 



k B T 
q-1 ' 



1, JSl yields the usual Boltzmann 



(7) 
(8) 



i.e., E can take any possible value if q < 1, but the distribution function has an energy cutoff if q > 1. Notice that 
both Q and hold for the extensive limit q = 1, in which case they give the expected constraint < E < oo. 

Due to the cutoff (J7J, if q > 1 f g is always normalizable. However, normalizability of the distribution function for 
< E < oo when q < 1 is not guaranteed. The integral 



f q (p)dp 



(9) 



is finite if f q (p) ~ p s when p — * oo, with s < — 1. Using (JSJ), 

hip) 1 



Cp'- 



where C = 



9-1 
k B T 



> 0. Thus, f q (p) is normalizable if q > 0. 



Let us now analyze the case q > 1. The cutoff in and © give 




(10) 



If 



ksT/ (q — 1) > mc , 



(11) 



then 111 Oil yields the cutoff in momentum space corresponding to the energy cutoff (J2J. In the limit case fcgT /((?— 1) = 
mc 2 , the only momentum allowed is p — 0. If k B T /{q — 1) < mc 2 , then no real value of p satisfies condition Q, and 
no distribution function exists. We may interpret this saying that for a given temperature, only q values consistent 
with lfTT)l are possible: 



q< 1 



k B T 



(12) 



3 



Summarizing: 

Case 1: < q. 

Forbidden q. f q is not normalizable. 

Case 2: < q < 1. 

f q is normalizable. Energy can take values 



< E < oo 



Case 5: 1 < <? < 1 



mc* 



f q is normalizable. Energy can only take values 



< E < 



k B T 
q-1 



i.e. 



< \p\ < 



k B T 
(9-1) 



Case A: 1 



Forbidden q. No values of p consistent with energy cutoff. 

Two main differences arise with respect to the non-relativistic case -9], where allowed values are — 1 < q < oo. In 
the relativistic case constraints are much tighter, < q < 1 + fc^T/mc 2 , and depend on the ratio of thermal to rest 
energy. 

Considering these constraints, the relativistic distribution function can be plotted. The result is qualitatively similar 
to the nonrelativistic case: 
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FIG. 1: Relativistic dispersion function © (feT = m = c = 1). Full circles: q = 0.5. Solid line: q = 1. Open circles: q — 1.5. 



III. DISPERSION RELATION FOR LONGITUDINAL PLASMA OSCILLATIONS 



We consider a proton-electron relativistic plasma. Assuming protons remain at rest, and that particles only move 
in one dimension (z), the dispersion relation for longitudinal oscillations is 



, ■< f, dfg(p z )/d Pz 

1 = i-f-m / d Pz — H — 

k J s + ikv z 



(13) 
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where lo p is the plasma frequency, k is the wavenumber, s is the argument of the Laplace transform (s = —iu, where u> 
is the wave frequency, if there is no damping), m the electron mass, p — rwyv is the particle momentum, and f q (p, z) 
is the equilibrium electron distribution function, normalized to unity. Particle energy is E = yfm 2 c A + pic 2 . 

Following Ref. 0, we now consider the long wavelength limit k <C kjj, where fco = 2n/X]j, and Ad is the Debye 
length. The integrand has a pole at v = —is/k. If Im t)<l, integration can be performed on the real axis: 

where p max is the maximum momentum allowed by (|10Jl if 1 < q < 1 + fcsT/rnc 2 , or oo if < q < 1. Since f q is an 
even function of p z , the first term in (|14|1 vanishes. The dispersion relation can then be written: 

muj l r max j. ^ 5 /? 



1 = — 2~ / dp* „ 

where 7 = (1 + pl/rn^c 2 ) 1 ^ 2 . Integrating by parts, and using the fact that / q (±p max ) = 0, finally yields 



Equation (|15(l depends on g in general, and therefore, unlike the nonrelativistic case [[j, nonextensive corrections 
to the dispersion relation of longitudinal waves may appear to the lowest order in kv z /s. 



IV. ULTRARELATIVISTIC PLASMA 

Closed expressions in terms of known functions can be obtained for Eq. (|15|l in the ultrarelativistic case \p z | ^> inc. 
As shown in Sec.^ q can only take values between and 1 + ^T/mc 2 . We consider two situations: 

A. < q < 1 

Since p max = 00, in this case, the dispersion relation Ijl5(l can be written 

ui f 00 1 1 
1 = - 2 ^./ ^372- (16) 



Noting that — (q — 1) > and — l/(q — 1) > 0, normalization of the distribution function 



UiPz) = A q — ^r: . (17) 

\-iq-Dm~*- 1 



k B T) 

I = I f g (Pz)d Pz , (18) 







yields 



where 



A q =— f , (19) 
mc 2r 



k R T 



The integral in l|16(l can be done analytically, yielding 



(20) 



1 = --£-,/ r, > ( 21 ) 

5 t V g — 1/ 



where 



J(r,a) 



-.{-Ta) a+1 aT 



V 7T 



1 — a 



3 + a 
2 

r 



r 

2 + a 



2^1 



1 + a 3 + a 3 



2 + a a 1 



' 2' 2' 



2 2 

-err 



1 — a 

2-Fi is the hypergeometric function, 



3F2 



2 2 



2 — a 3 — a 



-r 2 a 2 



2 F 1 (a,b,c;z) = ^ 



{a)u(b) v z v 
(c)„ v\ 



with [a) v = a(a + 1) • • • (a + v — 1), and m -F n is the generalized hypergeometric function, 

(oi)w • • ' (am)!/ 



jF n ({ai, . . . , a m }, {61, . . . , 6„}, z) = 



B. 1 < q < 1 + k B T/mc 2 
Now p max is given by Eq. ljTU|) , and dispersion relation Eq. lfH)|l is 



1 = - 2 Tf^ 



dp 2 



( 1 + i^f) 



3/2 



1 - (? - 1) 



Noting that q — 1 > 0, normalization of the distribution function yields 



2mc t 



i-(i-(g-l)J 



9 n -1 



where 



P 



Again, integration in (|25l) can be performed analytically, yielding 

, ,2 



1 



s 2 ^TTW 



where 



J(r,a) 



1-1 



a+l 



(a 2 r 2 + l) 2 



3 3 



ar ar 



(1 + iarr) y/T+0 2 ■ F x 1 + a; -, -; 2 + a 

2 2 ar + i ar — i 



1 



P 



l+a 



3 3 

[(l + mr)+p(i + aT)] • F x ( 1 + a; -, -; 2 + a 



where F\ is the Appell hypergeometric function of two variables: 

(&l)m(&a)n 



OO OO 



Fi(a;6i,6 2 ;c;x,?/) = ^ 



m— n— 



m! n! (c) 



a; m 3/ n 
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V. CONCLUSIONS 



Some consequences of describing a plasma based on a nonextensive statistical treatment have been discussed. In 
particular, the one dimensional plasma distribution function which maximizes Tsallis's entropy S q has been considered. 
In the nonrclativistic case, normalizability of the distribution function constrains q to values greater than or equal 
to —1. However, in the relativistic case a much shorter range is possible, and the maximum value of q depends on the 
ratio of temperature to particle rest energy. 

Then we derive the dispersion relation for longitudinal oscillations in a plasma composed of electrons and rest 
protons. In the long wavelength limit, and unlike the nonrclativistic case, nonextensive corrections appear to the 
lowest order. Finally, in the ultra-relativistic regime, this dispersion relation is written explicitly in terms of generalized 
hypergeometric functions. 
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